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Abstract 

The contribution of thermal fluctuations to the widths of isoscalar giant multipole resonances 
(GMR) in heated nuclei is studied. Starting from the collisional kinetic equation, it is shown that 
an additional contribution to the nuclear friction and the corresponding GMR widths arises due to 
the nonlinear dissipativity effect. It is also shown that the magnitude of the contributions of the 
thermal fluctuations to the nuclear friction coefficient and the GMR widths do not exceed ~ 20%. 
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I. INTRODUCTION 



In general, the damping of collective excitations in a cold finite Fermi system, e.g. the 
width of a giant multipole resonance (GMR) in cold nuclei, is determined by the two-body 
collisions, the particle-hole energy fragmentation (Landau damping) and the escape width. 
The damping in cold nuclei was intensively investigated in both the quantum (RPA like) 
and the semiclassical (kinetic theory) approaches □ □ □ □ □ 111 M2 . The 
collisional damping is due to the coupling of particle- hole excitations to more complicated 
states. In the kinetic theory, this type of damping is simulated by the collision integral and 
leads to a collisional component of the intrinsic width of the collective eigenstates. The 
fragmentation width is caused by the interaction of particles with the time- dependent self- 
consistent mean field. This contribution to the intrinsic width does not reflect a motion 
of the system toward the thermal equilibrium but indicates rather a redistribution of the 
particle-hole excitations in the vicinity of the collective state. In a hot system an additional 
contribution to the damping of the collective excitations arises through thermodynamic 
fluctuations of the corresponding collective variables because of the fluctuation-dissipation 
theorem. In this context, one of the most important open problem is the behavior of the 
GMR width in hot nuclei as a function of the temperature T. There are two essential 
different sources for the T-dependence of the GMR width. The first one is given by the 
thermal contribution to the damping width from an increasing nucleon-nucleon collision 
rate (2p2h excitations) plus a Landau spreading due to thermally allowed ph transitions 
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21, m 113 HJ, lia Ug . In the second one a temperature increasing of the width is caused by 



the coupling of the GDR to the thermal shape fluctuations of the nucleus 
the present work, we study a new effect of the influence of the thermal shape fluctuations 
of the nucleus on the damping of the collective motion caused by the nonlinear dissipativity 
appearing in the higher order variations of the collision integral. We point out that in the 
commonly used linear order of the variation of the collision integral, the thermal fluctuations 
do not lead to dissipation (viscosity) in the macroscopic equations of motion because the 
following ensemble smearing of the kinetic equation washes out the fluctuation terms from 
the final macroscopic equations of motion. This paper is organized as follows. In section II we 
suggest a proof of the Langevin equation for nuclear local variables (particle density, velocity 
field and pressure tensor), starting from the collisional kinetic equation. We perform a high 
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order expansion of the collisional integral and derive the non-Markovian pressure tensor 
to the Navier-Stokes-like equation of motion. In section III we carry out the ensemble 
averaging and reduce the local equations of motion to the macroscopic form and derive 
the macroscopic response function. In the derivations, the main features of the dynamical 
distortion of the Fermi surface are taken into account. Results of numerical calculations are 
presented in section IV. We conclude and summarize in section V. The Appendixes provide 
a derivation of the high order variations of the collision integral with respect to the variation 
of the phase space distribution function. 



II. LOCAL EQUATIONS OF MOTION 



We will start from the collisional kinetic equation in the following form 

df pdf dVdf - df rn rjP . 
at m or or op op 



(1) 



Here, / = f(r,p,t) is the Wigner distribution function, V = V(r,p,t) is the self-consistent 
mean field, _F ex t is the external driving force and 5St[/] is the collision integral which takes 
into account the memory effects 2l|. The random variable y = y(f, p, t) in Eq. (Q) represents 
the random force. As such, its ensemble average vanishes, (y) = 0. To reduce Eq. to 
closed equations of motion for the macroscopic collective variables we will follow the nuclear 



22 



23, 



24. 



25 



261 ] and take into account the dynamic Fermi-surface 



fluid dynamic approach 
distortion up to multipolarity 1 = 2. Evaluating the first three moments of Eq. in the p 
space, we reduce Eq. ((H) to the hydrodynamic-like equations of motion for particle density 
p, velocity field u and the pressure tensor P a/3 (for details, see Ref. j^ ) 

dp d 
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Here. 
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(p a -mu a )(pf3-mup)f, (5) 
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g = 4 is the spin-isospin degeneracy factor, Q Q/ g is associated with dissipative processes 

® af3 = ml (2nh) 3 <J)a ~ mUa ^ Pl3 ~ mu p) SSt \fh ( 6 ) 
and y a p gives the contribution from the random force 

1 r gdp 

In Eqs. (0) to (HJ) and in the following expressions, repeated Greek indices are to be under- 
stood as summed over. The pressure tensor P Q/ g can be written as 

Pa/3 — Pecfiafi + P'afii (8) 

where 

P ' aP = ml (2nh) 3 ~ mUa ^ p P ~ mu P) 6 hi ( 9 ) 
5f 2 = Sf 2 (r,p;t) is the small quadrupole deviation of the distribution function / from the 
one in equilibrium, / eq = f eq (r,p), and P eq is the pressure due to the Fermi motion of the 
nucleons in the ground state of the nucleus. Assuming the Thomas-Fermi approximation 
for /eq, one has 

eq 3 J (2nhf 2m /eq 5m{ g J Peq ' 
where p eq = p C q(r) is the particle density in equilibrium. Using Eq. (jHJ) we will rewrite Eq. 
flU) as 

d 

—P' ap + P eq A Q/3 + LP' a(i = Q af3 + y a/3 , (10) 



where 



a>3 drfs dr a 3 a/3 dr v ' 



and the operator L is defined by 



LP' -v 9 P' +P> du » +P> du P+ P > du « 



To evaluate t 
general form 



re tensor Q a p in Eq. (fTTH) we will use the collision integral in the following 



§m = J g^d^dp, wm)Qm) 5{Ae)6im (13) 

where w({pj}) = w(pi,p 2 ;P3,P4) is the spin-isospin averaged probability for two-body scat- 
tering, Q({fj}) = /i/ 2 (l - / 3 )(1 - / 4 ) - (1 - /i)(l - / 2 )/ 3 / 4 is the Pauli blocking factor, 



Ap = pi + f>2 — P3 — P4 and Ae = e\ + e 2 — £3 — £4, with ej = p|/2m + V(rj-) being the 
single-particle energy. In performing the variation of the Pauli blocking factor Q{{fj}) with 
respect to 5f , we will keep all the terms, up to the third order in Sf. The collision integral 
then takes the following form (see Appendix A) 

5St = <?Sti + 5St 2 + 5St 3 , (14) 

where SSt n is the variation of the collision integral SSt in the n-th order of 5f. Considering 
Eqs. (P) and (JHJ), the tensor Q a p can be written as 

where Q^) is due to the contribution of 5St n in Eq. (jBJ). The first order term Q^l is 
simplified as j3] 



q2 = (is) 



'a/3 _ J a/3' 

where r 2 is the two-body relaxation time in the case of quadrupole deformation of the Fermi 
surface. The higher order terms Q a L and Q a l can be reduced to the following forms (see 
Appendix A, Eqs. (1AT81) and (jCTjl ) 



n (2)_mP', ( 3 ) _ m 2 P> 2 

Va/3 _ 7 r aj3i V a (3 ~ ^ - r a/3; l 1 '/ 1 



where the quantities ( and £ are deduced from the collision integral (see Appendix A, Eqs. 
ifftTfljl and (533} ) and 



- ^(P'xx + P'yy ~ P'zz) ■ ( 18 ) 
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To simplify the derivations, we will introduce the operator as 



Equation (fTUj) is then rewritten as 

+ P eq A Q/3 + + ZP^ = A>P^ + y a/3 . (20) 

We will look for a solution of Eq. (j2*U|) in the following form 

p> _ p'(0) 1 p'(NL) / 91 \ 

^0/3 — ^a/? + ^a/? ■ l Zi J 
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iere, the tensor P'^g is obtained as a solution to the following linear differential equation 



f 



d p'(0) , p A , 1 p /(o) _ 

and it is given by the following non-Markovian form 

't'-t 



P$? (*) = -[* dt> expf- 

«/ — oo V 
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[P eq A^(t') - MO] • 



(22) 



(23) 



The tensor P^^ in Eq. (|21j) satisfies the nonlinear differential equation 



d_ p /(NL) 



l(NL) 

a/3 



LP' ap - NPL« = 0. 



a;3 



(24) 



dt af3 r 2 

To solve Eq. (J2H), we will use the iteration procedure. The first order iteration P^Bi it') to 
Eq. H2D reads 



dt' exp 



r 2 



(25) 



and the second iteration is given by 



dt exp 



(26) 



Below, we will apply Eqs. (J2J) to (jlj) to the small-amplitude vibrations of the particle density 
Sp near the equilibrium. We point out that we do not assume the velocity field u to be small. 
Finally, taking into account the above mentioned derivations we will rewrite Eq. (j5J) as 

5 2 e 



du r 



3u r 







Sp 2 



eq 



f^ + ^^-p eq iW = (27) 



<9r n 



<9r n 



where e is the particle energy density. 



III. ENSEMBLE AVERAGING AND MACROSCOPIC RESPONSE 

Let us introduce the displacement field x related to the velocity field u by w(r, t) = x(r, t) , 
where the dot denotes a time derivative. For the displacement field we will assume the 
following separable form x(r,t) = (3{t)v{f). Using this separable form of x(r, t), we reduce 
Eq. (j2*Tj) to the equation of motion for the collective variable j3{t) in the presence of the 
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external field J-' C y±{t) and the random force y(t) (see Appendix C, Eq. (|C1[1 ). Below, we 
will look for the response of a nucleus to the periodic external field J r e xt{t) = exp(iut). 
Because of the random force y(t) in Eq. (jClj) . we will separate the description of the 
collective motion into two parts with f3(t) = j3(t) + 5/3(t). The first motion is related to the 
driving force JF ext (t) and it is associated with the velocity (3. The second one is due to the 
random force y(t) with the velocity 5(3. We will assume that \5(3\ 3> \(3\. Performing the 
ensemble averaging, one can write 



P(h)P(t 2 )P(t 3 )) « P(h) {5(3(t 2 ) 5p(t 3 ))+/3(t 2 ) {SPih) 5p(t 3 ) / 

+ hs)(sp(h)SP(t 1 )). (28) 

We will also assume the following ergodic property for the correlation function 
(6p(t) 6p(f)), see Ref. 0, Ch.12, 

(Sp{t) 5$(t')) = r %L <y) e-^-*'\ (29) 

The Fourier component, (6f3 2 ^j , of the correlation function is governed by the correlation 
properties of the random force y(t), see below. The macroscopic equation of motion (jClj) 
is significantly simplified in the case of a Fermi distribution for the equilibrium distribution 
function f eq 



eq 



1 + exp 



t-t F 
T 



(30) 



where is the Fermi energy. In this case, one obtains from Eqs. (|A19|) that 1 (see 

also Figs. 1 and 2) and the contribution of the terms with A 2 , A 3 and A 5 in Eq. (jCl|) is 
negligible. Performing the ensemble averaging of Eq. (jClj) . using Eq. ()28j) and (y(t)) = 
and assuming (Pit)) = P(t) = p w exp(iut) , we reduce Eq. (|C1J1 to the following form 

AJ ( S P 2 ) 4t 3 Ai r°° du' ($P 2 ) 



1 + iur 2 \ J-oo 2tt 1 + iu'r 2 1 + iuT 2 

du' {Sp). 1 rlA, r du' {^) 




+ t 2 2 A 4 I — - '-f — + . 

-oo 2n 1 + iu'r 2 1 + i(uj' + uj)t 2 1 + iujt 2 J-oo 2n 1 + + uj)t 2 

(31) 

Considering the nuclear isoscalar quadrupole mode, we will assume an irrotational motion 
with the displacement field v(r) given by [31] 

v(r) =V (r 2 F 20 (f))/2, (32) 



and the time dependent radius of the nucleus given by 



R(t) = R l+p(t)Y w (r) . (33) 



In this particular case the calculation of the coefficients A\ and from Eqs. (jC4j) and (jC7|) 
gives: A 1 = 16 A (mP cq ) 2 / £ and A 4 = 12A . The mass coefficient B of Eq. (JOS) for the 
displacement field of Eq. (f3l2j) is given by 

B = ^-AmRl (34) 

where A is the nuclear mass number. Let us introduce the collective response function x( a; ) 
as 

= Tw (35) 
Using Eqs. (}3~T]) and we obtain from Eq. (|3l)Jl 

X -\uj) = -u 2 + cu 2 + iu 7o + 12r 2 2 ^ ^ AT(^). (36) 

D 1 + ILUT2 

Here, 



and 



" = V 5 ' 7o= 5-TT^F (37) 

7 -00 27r 1 + zcj't 2 I 1 + z(u/ + a;)r2 1 + iur 2 1 — ^'r 2 

+ — 1 Vf T \ ) • (38) 

1 + iu)T 2 1 + i(u' + uj)t 2 J 

We have also used the following notations 

C '^ = A °TT& (39) 

and 

t" = £/(mP cq ) 2 . 

We point out that the additional contribution to the stiffness coefficient C'(ui) in Eq. (071) 

is caused by the distortion of the Fermi surface |29|. The expression (jHBj) can be rewritten 

as 

X -\uo) = (oj 2 + Alu 2 - lu 2 ) + iu ( 7o + A 7 ), (40) 



S 



where we have introduced the following notations for the additional components of the 
relaxation coefficient and the squared frequency uj 

— = 12r 2 2 {ReK(u) + uT 2 lmK(uj)} , (41) 

7o 



f = &,{B*KM-±l^ lu) }(l + °ffi)-\ (42) 



Note that above, u is real. Finally, the macroscopic strength function S(u) = -Imx(w) is 
given by 



S ( U ) ' .? , a. ,2 _ ,.,2^2 i i A^2,.,2' ( 43 ) 



(70 + Aq/)u 
(oo 2 + Auj 2 - uo 2 ) 2 + ( 7o + A 7 ) 2 ^ 2 
Both the additional spreading A7 and the resonance shift Aujq appear in the strength 

function (}4*3"j) due to the nonlinear dissipativity effect. 



IV. NUMERICAL RESULTS AND DISCUSSION 



We have performed the numerical calculations assuming a Fermi distribution for the equi- 
librium distribution function of Eq. (J3U)) and adopting the Fermi energy ep = 39 MeV and 
the nuclear radius Rq = tqA 1 ^ 3 with tq = 1.12 fm. The higher order relaxation parameters 
t' = (/mP eq and t" = £/(mP eq ) 2 are related to the collision integral and can be evaluated 
using Eqs. (jA19|) and (jA31|) from Appendix A. We point out that in the limit of a cold 
nucleus, T — > and / eq = 0(ef — e), the corrections r' and r" take the following simple 
form 

1 1 m 3 P 2 w , . 

0, f-^, 44 

where pp is the nucleon Fermi momentum and the scattering probability wq = 15ir 2 h 5 /m 3 ga 



is related to t 
MeV from 



re in-medium cross section a- m of nucleon-nucleon scattering. We use a = 9.2 
which corresponds to cr in (Tf rcc /2, where <Tf ree ~ 40 mb is the cross section 
for the nucleon-nucleon scattering in free space. Note that both relaxation parameters r' 
and t" can not be directly interpreted as the corrections to the observable relaxation time. 

In particular, the value of 1/ r" does not equal to zero in the ground state of the nucleus. The 

(2) 

relaxation parameters r and r determine the contribution of the viscous tensors Q a L and 
Qal ( see Eq- (fTTjl) to the local equations of motion and both tensors Q^p and Q^l disappear 
in the ground state. We also point out that the relaxation parameters r' and r" as well as 
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FIG. 1: Temperature dependence of H/t' (solid curve 1) and h/r" (solid curve 2) for the case of 
the temperature-dependent Fermi distribution function ()^()|) . The dashed line is the calculation of 
h/r" with the sharp Thomas- Fermi distribution function 0(e^ — e). 

( and £ depend on the nuclear mean field potential V due to the space integrals r^- and 
Tijk, see Eqs. (|A10|) and (jA21|) . This dependence appears after the Abrikosov-Khalatnikov 
transformation (jA7|) in the collision integral 5St[/]. However, due to the presence of the 
strongly picked functions df^i/dei, at e = €p, in Eqs. (|A10|) and ()A21|) the final results for 
( and £ are not sensitive to the specific choice of the mean field potential V at T <C ef- In 
Fig. 1, we have plotted the results of calculations of the quantities h/r' (solid curve 1) and 
h/r" (solid curve 2) as functions of temperature, T, for the nucleus with A = 224. 

Both quantities h/r' and H/r" show a very broad and weak maximum. The magnitude 
of the maximum does not exceed the value of 0.07 MeV for h/r' and 0.49 MeV for h/r". 
The dashed line in Fig. 1 corresponds to the value of h/r" from Eq. (J44j) . We can see from 
Fig. 1 and Eq. (|44|) that the simplest Thomas- Fermi distribution function 0(e^ — e), with 
H/r' and h/r" from Eq. (|44|) . provides a good description of both quantities h/r' and h/r". 
Fig. 2 shows the ratio of the collisional relaxation time Ti to both relaxation parameters r' 
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FIG. 2: Same as in Fig. 1 but for the ratio tz/t 1 multiplied by the factor 10 (solid curve 1) and 
the ratio Tzjr" (dashed and solid curves 2). 

and r" . For the relaxation time r 2 , we have used the expression from Ref. ^| which takes 
into account the memory effects. Namely, 



r-2 



(45) 



(^ ) 2 + 4tt 2 T 2 ' 

As seen from Fig. 2, the value of tz/t' is relatively small over the entire range of the 
temperature. The value of t 2 /t" decreases with the temperature monotonically starting 
from T2/t" = 1.16 at zero temperature. 

Let us now carry out a numerical study of the additional contribution to the friction 
coefficient, A7, caused by the nonlinear dissipativity, see Eqs. (J4Ti|) and (jH)) and the corre- 
sponding contribution to the width T of the isoscalar giant quadrupole resonance (GQR). 
To apply Eqs. (|3*Hj) and ([41)1. we have to derive the spectral correlation function {8(5' 2 
Using the correlation properties of the random force [30] 

27oT 



B 



(46) 
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we obtain according to the fluctuation-dissipation theorem the following result 

2Vu 2 
w (ul - uj 2 ) 2 + >y$uj 



where Z> is the diffusion coefficient 



2? = ^. (48) 



To evaluate the relative contribution to the collisional width T from A7, we will start from 
the usual case with A7 = 0. In this case, the width T of the GQR can be obtained from the 
solution in the form u) = Retu + iT/2h to the following secular equation, see Eq. (J43|) . 

K 2 -c 2 ) 2 + 7o V> = 0. (49) 

For the numerical solution of Eq. (jjQj) , we have used in Eq. (pTTj) the liquid drop stiffness 
coefficient Cldm hi the form |33| 

c LDM = _L (i _ 1)(£ + 2)M * _ AAzi 6c |l, (50) 

where b s = 17.2 MeV and be = 0.7 MeV are respectively, the surface and Coulomb energy 
coefficients appearing in the nuclear mass formula. 

Fig. 3 shows the results of the numerical solution of Eq. ()49|) for the nucleus with 
A = 224 and Z = A(l - 6 • 10~ 3 A 2/3 )/2 = 87, which corresponds to the valley of beta- 
stability 33]. The energy of the collective excitation E = h Re u decreases with temperature 



and approaches the hydrodynamic (liquid drop model) limit Eldm = Ti\J Cldm / B at high 
temperatures. In Fig. 4 we have plotted the temperature dependence of the parameter 
ET2/T1 which determines the sound regime: ET2/T1 ^> 1 for the zero sound (rare collision 
regime) and ET2/h 1 for the first sound (frequent collision regime). 

The solid curve in Fig. 4 corresponds to the calculation with the temperature dependence 
of E given by Fig. 3. For the dashed line, the phenomenological parametrization of the GQR 
energy E = E R = 60 • A~ 1/3 MeV was used. Using Eqs. (jSHJ), (J4~T|) . (|4"7j) and ()48j1 . one can 
evaluate the contribution A7 to the friction coefficient due to the nonlinear dissipativity. 
In Fig. 5, the value of A7/70 is shown as a function of temperature. The ratio A7/70 
equals to zero at T = because A7 appears due to the thermodynamical fluctuations of the 
collective variable (3. The ratio A7/70 increases with temperature and reaches a maximum 
value, which does not exceed ~ 0.2. 
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FIG. 3: Energy E = /iRew of the isoscalar giant quadrupole resonance (GQR) and the corre- 
sponding collisional width V for the nucleus with A = 224 as obtained from Eq. (|49[) . 

In the high-temperature region the ratio A7/70 decreases because the temperature de- 
pendence of 70 ~ T 2 is stronger than that of A7 ~ T. For comparison, we have also 
performed the calculation of the ratio A7/70 using the phenomenological parametrization 
for the GQR energy Er = 60 • A^ 1 ^ 3 MeV (see dashed line in Fig. 5). In this case, the vari- 
ation of A7/70 with temperature is somewhat stronger. Taking into account the nonlinear 
dissipativity effects, the collisional width V of the GQR is obtained from the solution, in 
the form uj = Reu + iT' /2h, to the secular equation, see Eq. (J4"H|) . 

(cu 2 + Alu 2 - oo 2 ) 2 + (70 + A 7 ) V = 0. (51) 

In Fig. 6 we have plotted the temperature dependence of the widths T (dashed lines) and 
T' (solid lines) for two choices of the resonance energy: E = hKeu using Eq. (J49)) (curves 
1) and E R = 60 ■ A~ 1/3 MeV (curves 2). 

We point out that an increase of the width is more apparent for curves 1 in Fig. 6 because 
of the temperature dependence of E. The comparison of the solid and dashed lines in Fig. 
6 shows that the contribution of the nonlinear dissipative effects to the width T does not 
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FIG. 4: Dependence of the dimensionless parameter Er2/h on the temperature T for the GQR in 
the nucleus with A = 224 with T2 from Eq. (|45j). The solid curve was obtained using E = ?iReu; 
from Eq. (|49|h The dashed curve was obtained with E = En, where Er = 60 • yl _1//3 MeV is the 
experimental value of the GQR energy. 

exceed ~ 20%. 

In Figs. 7 and 8 we have plotted the strength function S(u) of Eq. (|4"3j) . The comparison 
between the solid and the dashed lines in Fig. 7 shows the accuracy of the derivation of the 
value of A7 directly from the strength function S(oj) of Eq. (J43j) and through the solution 
of the secular equation (|51|). The comparison of the solid lines with the dashed lines in Fig. 
8 demonstrates the effect of the nonlinear dissipativity on the strength function. 



V. SUMMARY AND CONCLUSIONS 



Starting from the collisional kinetic equation with a random force and using the p- 
moments techniques, we have derived the equations of motion of the viscous fluid dynamic 
for the local values of particle density, velocity field and pressure tensor. The obtained 
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FIG. 5: Temperature dependence of the ratio A7/70 for the nucleus with A = 224 for the GQR. 
The solid curve was obtained using Eq. (|41|) with uj from Eq. (|49[) (see also Fig. 3); the dashed line 
was obtained using Eq. (|41j) with 70 from Eq. ()37|) and u = ajR = Er/H with the phenomenological 
parametrization Er = 60 • A^ 1 ^ MeV. 

equations are closed due to the restriction imposed on the multipolarity / of the Fermi sur- 
face distortion, up to Z = 2. The important features of these equations of motion are due 
to the non-Markovian form of the pressure tensor P a p. In contrast to the commonly used 
r-approximation, we take into account the higher orders of the variation of the collision 
integral with respect to the variation of the phase-space distribution function. Using the 
Abrikosov-Khalatnikov transformation we have then obtained the collision integral in the 
form of the extended r-approximation. Assuming a separable form for the displacement 
field, we have introduced the macroscopic collective variable f3(t) and reduced the problem 
to a macroscopic equation of motion for (3{t). Note that we do not assume the velocity (3(t) 
to be small. The final macroscopic equation of motion (jCl|) includes both the memory ef- 
fects and the nonlinear dissipativity terms ~ f3 3 . We have separated the description of the 
collective motion into two parts. The first (slow) one is related to the driving force Fextit) 
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FIG. 6: Collisional width T as a function of temperature T for the nucleus with A = 224 for 
the GQR. The solid lines are for V = 2fiIm.Lo from Eq. (|49j) and the dashed lines are for V 
from Eq. (j51|) . The curves 1 were obtained using the temperature dependent resonance frequency 
u = ojr = Re a; from Eq. ()49|) . The curves 2 were obtained using u = ujr = Er/1% with the 
phenomenological parametrization En = 60 • A -1 / 3 MeV (see also Fig. 5). 

and it is associated with a slow motion having the velocity (3. The second (fast) one is 
due to the random force y(t) with the velocity 5(3 3> (3. Using the correlation properties 
of the random force, we have performed the averaging of the macroscopic equation of mo- 
tion over the fast fluctuations ~ 5(3, reducing the nonlinear dissipativity terms to the form 
~ (3{t) (5(3{t')5f3{t")^ , which is linear with respect to the slow collective motion ~ (3. Fi- 
nally, assuming a periodic driving force .Fext(*) ~ exp(icut), we have derived the macroscopic 
strength function S(u>). As seen from Eq. the nonlinear dissipativity effect leads to the 
additional spreading A7 and the resonance shift Ao; in the strength function S(u>). The 
contribution A7 appears due to the thermodynamical fluctuations of the collective variable 
f3. In contrast to the Fermi-liquid friction parameter 70 with 70 ~ T 2 (at T C ep), the 
spreading A7 is a linear function of the temperature T. This fact provides a non-monotonic 
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FIG. 7: The strength function S(uj) in ?i 2 -units for two temperatures: T = 1.9 MeV (curves 1) and 
T = 5 MeV (curves 2). The solid curves 1 and 2 were obtained from Eqs. (|41|h (|42|) and (|4H[) . The 
dashed curves 1 and 2 were obtained using Eq. (05J) with 70 + A7 = T/n and + Aw 2 , = {E/Kf, 
where E and T are obtained from the solution, in the form u = E/h + iT/2h, to the secular 
equation (|5Tj) . 

behavior of the ratio A7/70, see Fig. 5. As seen from Fig. 5, the nonlinear dissipativity 
effects are enhanced at the moderate temperatures T « 2 MeV and do not exceed ~ 20%. 
The nonlinear dissipativity effect increases the collisional width of the GMR. Usually the 
total collisional width of the isoscalar GQR in cold nuclei does not exceed 30-40% of the 
experimental value and the main contribution to the width is due to the Landau damping. 
One can expect that the nonlinear dissipativity effect on the collisional width can lead to a 
deviation of the temperature dependent width T(T) from the usual Fermi liquid prediction 
T(T) ~ T 2 . Unfortunately, at present time, experimental data on the temperature behavior 
of r(T) of the isoscalar GQR are not available. In this respect, it is more instructive to 
study the isovector giant dipole resonance where the temperature dependence of T(T) was 

. However our final results for the viscous tensor Q a p 




studied for some heavy nuclei 36, 
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FIG. 8: The strength function S{oj) in ?i 2 -units for two temperatures: T = 1.9 MeV (curves 1) and 
T = 5 MeV (curves 2). The solid curves 1 and 2 are the same as in Fig. 7. The dashed curves 1 
and 2 were obtained from Eqs. I|43|). but with A7 = 0. 

and the relaxation parameters ( and £ can not be applied directly to the isovector mode 
because the dipole distortion of the Fermi surface must be taken into account in the collision 
integral (|13j) . in contrast to our case of the isoscalar GMR, see Sect. II. The generalization 
of our approach to the case isovector modes is now in progress. 
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APPENDIX A 



As a basic expression for the collision integral 5St[/] we use Eq. (fT3j) . The second and 
third variations of Eq. (fT3*j) with respect to 5f take the following form 

g 1 ] dp 2 dp 3 dp<i 8 2 Q 



g 2 dp 2 dp 3 dp i ^ 5 3 Q 



6f(i)6f(j) 5(Ae)5(Ap), (Al) 



cq 



5st 3 = y ^(fe»E 



5f(t)5f(j)5f(k) 5(Ae)5(Ap), (A2) 

eq 



(2^)6 vu-jj^ (j)5/(A;) 

where = 5f(pi) and the symbol J] means a summation over indices i, j, A; = 1 -r-4 with 
i / j, j / I; ^ i. We will follow the fluid dynamic approach and represent the variation 
of the distribution function 5/ in the following form: 

df ■ 1=2 
Sf(i) = "^rMO, "(0 = E^ m; (r,t)F 2mi (a). (A3) 

We point out that the / = and 1 components of the expansion (|A3|) do not contribute to 
the collision integral (|13|) . reflecting the conservation of particle number and momentum in 
a collision. The expansion coefficients z/ 2m (r, t) in Eq. (jA3|) are related to the pressure tensor 
P' a p of Eq. ©. Using Eqs. @ and (jUJ), we obtain 

5 / 2 



where 



dpp 4 ^, (A5) 

and p = /p is the unit vector. In particular, performing the angle integration in Eq. (|A4|) . 
we obtain 

"-"i^M- (A6) 

where Pq is given by Eq. (fIS)l. To evaluate the collision integral <5St 2 , we will substitute Eq. 
(IA3|) into ()A1|) and make use of the Abrikosov-Khalatnikov transformation in the following 
form U 

J dp 2 dp 3 dpi (...) 5(Ap) ^~yJ v de 2 de 3 de 4 J (•••) > ( A7 ) 

where <if2 = sm8d6d(p, 9 is the angle between pi and p2 ; is the angle between the planes 
formed by {pi,p 2 ) and (p3,Pa), and <p 2 is the azimuthal angle of the momentum p 2 in the 
co-ordinate system with z-axes along p\. We point out that the angle <fi varies only from 
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to it because the particles are indistinguishable. Using the transformation (jA7|) and the 
relation (see Appendix B) 

r 2ir 



we obtain 



f Z7T d(p2 

/ -^—Y nm (fti)Y n/rn/ (Qj) = Y nm (Q 1 )Y n/m/ (Q 1 )P n (cos9 i )P n/ (cos9 j ), (A8) 

JO 271 



(27i) 2 m 3 — 
6St2 = 9 ~^Kf Hl)]2 E(^(cos^)P 2 (cos^))r^, (A9) 



Here, is given by 

S 2 Q 



de 2 de 3 de 4 



dti den 

cq ' J 



5(Ae), (A10) 



and the symbol (...) denotes the following average 

(w(9,(l))P 2 (cos9 i )P 2 (cos9 j )) = 2 f d9sin(9/2) (* ^w(9, 0) P 2 ( cos 9i)P 2 ( cos 9 j), 

Jo Jo 2ix 

where cos0j = (pj ■ p t ), i.e. 9 2 = 9, and 

cos 9 3 = cos 2 (9/2) + sin 2 (0/2) cos0, cos# 4 = cos 2 (9/2) - sin 2 (9/2) cos0, (All) 
and P;(cos#) is a Legendre polynomial. Using Eqs. (jOJ), (Jill), (|THj) and (|A9J) . we obtain 

= ^ E( wP 2(cos g i )P 2 (cos | dfii p 1)0 Pi^Kl)] 2 , (A12) 

where 

POD 

Rij = / dpipjnj. 
Jo 

To exclude the unknown amplitude z/(l) from Eq. (jA12|) . we will calculate the arbitrary par- 
tial contribution to the tensor Q^. Using Eq. (|A1|) . we will consider the partial contribution 
Qq/3,12 to tne tensor Q^p given by 

Qi 2 j,i2 = m ^)9 / dPiPi.cxPhpdfi 1 ) J dp 2 5f(2) J dp 3 dp 4 



' r ' 5/(1)5/(2) 



<S(Ae)<S(Ap). (A13) 

eq 
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We will assume the isotropic probability scattering w(9, <fi) = wq, and apply the Abrikosov- 
Khalatnikov transformation in the following form 

/vr? r 

dp 3 dp 4 (...) 5(Ap) -— / de 3 de 4 # (...) , (A14) 

where pp is the Fermi momentum. Using Eqs. (JA3j) and (JA14j) . we transform Eq. (JA13J) as 

(2) g 3 27rm 3 w 2 - Y 2m {Vl 2 ) r . A 

g ^> 12 = m WW 12 Jk 2 2m l dn2 ^m)J dni Pl ' a Pl >" (A15) 

Integrating over dfl 2 in Eq. (|A15|) and using Eq. (|A6j) . we obtain 

(2) 2 37rm 3 w P i2 , f A A 

^' 12 = 5 (2?r^) 6 7 7 1 Pl,/3 ' ^ ^ 

Comparing Eq. (|A16J) with the partial i = 1, j = 2 term of Eq. (|A12|) . we find the relation 

/ dOi Ki)] 2 = 47rg( p 2 3 (cosg)) (27r / )3 wt ^ / ^ Pik. Pi^J Ki)- ( A17 ) 

Finally, from Eqs. flESJ), {23} and fATfj) we obtain 



Q$ = (M8) 

where 



3irgm 3 w E(^2 (cos #j)P 2 (cos 9j))Rij 



(A19) 



C (2tt^) 3 (P 2 (cos6))P 
Let us go now to the third order variation of the collision integral 5St 3 of Eq. (|A2[) . Applying 
Eqs. flA3|) and the transformation (jA7|) to Eq. (jA2j) and using the relation (see Appendix 
B) 

= F ri , m (ni)F n ' m '(ni)F„» m »(ni)P n (cos6'i)P n '( cos 8j ) P n » (cos 9k ) , 
we will reduce the collision integral <5St3 to the following form 

SSt 3 = -9 J h y E( wP 2(cos^)P 2 (cos%)P 2 (cos^))r^, (A20) 



where 

5 3 g 



de 2 de 3 de 4 



5f(i)5f(j)5f(k) 
Using Eqs. (ffij). ([TI)l . (fTHj) and (lA"20i we obtain 

1 f gdp 



dfcq,i dfeqj <9/ e q,fc x( « ^ /a<ii\ 



= — / (2^)3 ^° ~ mM «)(^/3 - ^M/3)5St 3 
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g 3 (27r) 2 m 3 ^ 



^2(wP 2 (co S e i )P 2 (cose j )P 2 (co8e k ))R ijk JdQ 1 p ha p 1 ^ [z/(l)] 3 , (A22) 



m (2vrft) 9 

where 

Rijk = I dpp 4 r ijk . (A23) 



oo 

4„ 



(2) 

Similar to the previous evaluation of the tensor Q a p, we will consider the partial term of 
Eq. ()A20|) with % = 1, j = 2, k = 3 and the corresponding partial tensor Q^l 123 which is 
given by 

x | dpW(2) | dpW(3) | dp 4 w(e,<l>)6(Ae)6(Afi, (A24) 
where we have used the following relation 

5 3 Q 



= -1. 

eq 



5f(l)5f(2)5f(3) 

We will again assume the isotropic scattering probability: w(9,(/>) = wo and apply the 
transformation (jA14J) to Eq. (jA24|) . The angle integrals over dQ 2 d<p, appearing in Eq. 
(1A24J) . can be transformed as 

/ dn 2 d / 2m '^ Y 2m// (n 3 ) = / dn 2 def> p ^ c ™^ Y 2m ,(n 2 ), (A25) 

J cos(^/2) J cos(6y2) 

where we have used the relation j3| 

r ^Ynmi^i) = F„ m (fii)P„(cos^). (A26) 

JO 2,71 



The result reads 



(3) g 3 m 3 w 2 

^,123 " --2(^)9^3^^ 



x I dn 2 d<p P ^^ Y 2m (n 2 ) J dfij p 1)Q p 1)/3 [z/(l)] 2 . (A27) 

Performing the integration over dfi 2 d0 in Eq. (gZfl and using Eqs. (|A17|) and (|A6Jl . we 
obtain 

n(3) 27 gm 3 w P 123 2 /• a . m , AOfiA 

^ = -28 (2^) 3 (P 2 (cosg)) ^ mP ° 7 d ^ Pl '° V ^ V[X) - (A28) 

Comparing Eq. (|A28|) with the partial % = 1, j = 2, k = 3 term of Eq. (|A22|) we obtain the 
following relation 



/ dfii pi )tl pi )(S z/ 3 (l) 



27 (2tt?1) 6 m 2 P ( 



2 c>'2 




28 (2vr)V^ 2 (P 2 (cosfl))(P 2 (cos#)P 2 (cosfl 3 )> 
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x / dQi p 1>a p ltP 1/(1). (A29) 
Finally, substituting Eq. (jA29j) into Eq. ()A22|) and using Eq. (J23J), we obtain 



Q$ = (A30) 



where 

1 = 27m 3 w E(P 2 (cos ^)P 2 (cos fl 3 -)P 2 (cos fl fc ))P, jfc 
e 28 (P 2 (cos#))(P 2 (cos#)P 2 (cosfl 3 )>/ 3 ' 



(A31) 



APPENDIX B 



In this Appendix, we will consider some angle integrals which appear in the calculations 
of the collision integral and its variations. Let us start from the integral 

M 23 = J ^Pr(cose 2 )e im * 2 P^'(cose 3 )e im * 3 P i (cos^), (Bl) 

where (Qj,<&j) are the angle coordinates of the momentum vectors pj in the arbitrary co- 
ordinate frame (j = 1 -j- 4) and 9 is the ang le between the vectors p\ and p 2 . Using the 
addition theorem for spherical harmonics [35] 

^(cosfl) = £ ^!l!pi r| (cose 1 )pi'-i(c se 2 ) e ^-^, (b 2 ) 

r=-l i}+\ r \V- 



we find 



M- 



23 — z 

r- 



(l- 


\r\)\ 


(1+ 


\r\)\ 



x P\ A (cos Q 2 )P\ A (cos Q^e^e^e^i-^) 

= -^-P„ m (cos0 1 )e im * 1 Py(cose 1 )e jm '* 1 P n ,(cos0 3 ). (B3) 
2n + 1 

Here, 6 3 is the angle between pi and p 3 . On the other hand, using the direction of p\ as 
polar axis with dVL 2 = sin 6d6d<f)2 where 2 is the azimuthal coordinate of p 2 in the new 
'-ordinate frame, we will rewrite Eq. (JB1)) as 

M 23 = J y sin^(cos^) ||^P n m (cose 2 )e m * 2 Py(cos0 3 )e im '* 3 | . (B4) 



a, 
co 
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Here and below the angles 0j and <3>i are dependent on the angles 6 and 2 . Comparing Eqs. 
()B3|) and (jB4|) and using the orthogonality condition for the Legendre polynomial Pi(cos8) 
in Eq. (jB4|) . one obtains the following integral relation 



/ ^ j P„ m (cose 2 )e fm ' I>2 Py(cose3)e im '* 3 

= P n m (cose 1 )e im * 1 Py(cose 1 )e im ' <I>1 P n (cos^P ri .(cos^)- (B5) 
Starting from the integral 

M 24 = / ^P n m (cose 2 )e 4m ' I>2 Py(cose4)e im * 4 P / (cos^, 

J 4:71 



we will also obtain an integral relation analogous to Eq. (jB5|) but with the replacement 
3^4. Let us consider now the integral 

M 34 = / ^P n m (cose 3 )e im<I>3 P n 7'(cose4)e im ' <I,4 P(cos^)- (B6) 
J Am 

Using the addition theorem for p(cos6 , 3 ) (see Eq. (|B2j) ). we reduce Eq. (|B6|) as 

M 34 = P„ m (cose 1 )e ^m<I,1 / ^P™'(cose 4 )e m '* 4 

In + 1 J 2ix 



= ^ T Pr(cose 1 )e im * 1 P^'(cose 1 )e im '* 1 P n (cose 4 ). (B7) 

Replacing in Eq. (jB6|) the integration over $ 3 to integration over $ 2 and using the direction 
of pi as a polar axis, we will rewrite Eq. (|B6|) as 

M 34 = J ^sin^i(cos^) ||^Pr(cose 3 )e^Py(cose4)e m '* 4 | . (B8) 



Comparing Eqs. (|B7jl and (|B8|I and using the orthogonality conditions for the Legendre 
polynomials, we obtain 

^P„ m (cose 3 )e im<I>3 P™'(cos 0±)e im '^ 



2tt 

= Pf(cose 1 )e im * 1 P n Y(cosei)e 4m '* 1 P„(cos^ 3 )Pn'(cos^). (B9) 
Using the representation of the spherical function Y nm (Q) via the Legendre polynomials 

p™(cosfl) \m 



. 2n+l(n - m)! os 
\ 4vr (n + m)! 
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and collecting Eqs. (jB5|) and (jB9|) we obtain the following integral relation 

/ ~^ Ynm ^ Yn ' m '^"> = Y n m (^l)y n ' m '(^l)Pn{cOS 6 i)P nl {COS 6 j), (Bll) 

where i,j = 1 -j-4. Let us consider finally the integral 



M 234 = / ^P™(cose 2 )e im * 2 P n Y(cose 3 )e m,I>3 P n T"(cose4)e im * 4 P i (cos^. (B12) 

J 4:71 

Similar to the previous consideration, we will transform Eq. (jB12|) as 

f d6 

M 234 = J ysin0Pj(cos0) 



X 



^P n m (cose 2 )e im<I>2 P^'(cose3)e im ' ,I>3 P^"(cos Q A )e im '^ 



2n 

-P^(cosei)e iw01 Py(cosei)e im ' 01 



5nl 



2n + l 

x P n T,"(cose 1 )e im "* 1 P n ,(cose 3 )P„»(cose 4 ). (B13) 
Using the orthogonality conditions for the Legendre polynomials, we obtain 

/ ^P n m (cose 2 )e ^m * 2 Py(cose3)e ^m ' <I,3 P^;'(cos04)e im '' <I,4 
J 2rr 

= P™(cosei)e im * 1 Py(cose 1 )e im '* 1 P n T,''(cose 1 )e im " ,I>1 

x P n (cos 9) P n , (cos 9 3 )P n „ (cos 6±). (B14) 
Finally, taking into account Eq. (jBlOj) we will generalize Eq. (jB14|) as 

"^rrnit^il^i'm' (S^'j) Y n"m" (^fc) 



2vr 

F„ m (fii)y„/ m /(f2i)y ri » m «(f2 1 )P n (cos^)P^(cos0 i )P„»(cos^). (B15) 



APPENDIX C 

In this Appendix we give a proof of the macroscopic equation of motion for the nuclear 
shape variable (3(i) derived by the displacement field as x(r, t) = j3(t)v(r), see Sect. III. 
Substituting this separable form in Eq. (|2T|) and multiplying by v a , summing over a, and 
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integrating over r space, we obtain the equation of motion for the collective variable f3(t). 
Namely, 



B(3 + D /3 2 + A dt'exp 



— oo J — oo 



Di I dt' I dt[ exp 



${f) + C LDM f3 



r-2 



D 2 / dt' / dt' x j dt' 2 exp 



't[ - t s 



r-2 



exp 



t' 2 -r 



T 2 



/t ft' ft' ft' ( 

di! \ dt[ / dt' 2 / dt' 3 exp 
-00 J— 00 J —00 J— 00 V 



t[ - ? 



x exp 



't' 3 -t r 



T 2 

$&)$(&)$(&) 



exp 



. T 2 , 



+A 2 I dt' I" dt[ I 1 dt' 2 I dt' 3 exp 

00 



- r 2 , 



exp 



» T 2 . 



+A 3 / dt' I dt[ I 1 dt' 2 I 1 dt' 3 I dt^exp 

00 



. T 2 . 



exp 



t' —t'\ 



T 2 



x exp 



- T 2 , 



$(t' 2 )$(t' 3 )$(t' 4 



+A 4 / dt' / / d4 ex P 



t' 2 -f 

. T 2 , 



p{t')KAW 2 , 



+A 5 I dt' I dt[ I 1 dt' 2 l h dt' 3 exp 



% - V 

. T 2 , 



exp 



't' — t' " 

L 3 t 1 



T 2 



= BF cxt (t) + By(t), (CI) 

where BJ-" ext and B y(t) are, respectively, the external and random forces in the collective 
space of the variable (3 (we have separated the mass coefficient B from the external and 
random forces for technical convenience). The transport coefficients in Eq. (jCl|) are given 
by 



B = m J dfp eq v 2 , C LDM = J dr (J^j 



eq 



9 ( 



1 2 



Dn = m 



J drp eq v a i 



dr,. 



Pi 



dr -^A n A. 



D 1 = J dfi(P eq A, 
dv a 



dr,. 



--Ql\. au 

C or v 



A 



— J drP cq A au — 



Pi? 



A 1 nr j dr -^-A 2 Q A au — , 



(C2) 

(C3) 
(C4) 
(C5) 
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with 



and 



P /— — — — — — \ dv 

A 2 =m I dr-^ (A Qi ,L(P eq A ) + A L(P eq A Qi ,) 



Pl^ dv n 



A* = 2m 2 / df 



dr v ' 
j _x dv. 



A 4 = J drL (P cq A 



C 2 av dr v ' 
A 5 = m J dft( ^A aJ ^. 



L ^ " v " 9r, + ^ dr v + ^ 9r, + ^ 9r, ' 
x _dv a dvp 2 du„ 
dr^ <9r Q 3 0r„ 



Aq — ~{^xx + A ra — A zz ). 



dr v ' 



(C6) 
(C7) 
(C8) 
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